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In the standard approaches to neutrino transport in the simulation of core-collapse supernovas, one will often 
start from the classical Boltzmann equation for the neutrino's spatial, temporal, and spectral evolution. For each 
neutrino species, and its anti-particle, the classical density in phase space, or the associated specific intensity, 
will be calculated as a function of time. The neutrino radiation is coupled to matter by source and sink terms on 
the "right-hand-side" of the transport equation and together with the equations of hydrodynamics this set of cou- 
pled partial differential equations for classical densities describes, in principle, the evolution of core collapse and 
explosion. However, with the possibility of neutrino oscillations between species, a purely quantum-physical 
effect, how to generalize this set of Boltzmann equations for classical quantities to reflect oscillation physics 
has not been clear. To date, the formalisms developed have retained the character of quantum operator physics 
involving complex quantities and have not been suitable for easy incorporation into standard supernova codes. 
In this paper, we derive generalized Boltzmann equations for quasi-classical, real-valued phase-space densities 
that retain all the standard oscillation phenomenology, including the MSW effect (matter-enhanced resonant 
flavor conversion), neutrino self-interactions, and the interplay between decohering matter coupling and flavor 
oscillations. With this formalism, any code(s) that can now handle the solution of the classical Boltzmann or 
transport equation can easily be generalized to include neutrino oscillations in a quantum-physically consistent 
fashion. 

PACS numbers: 14.60.Pq, 05.30.Ch, 05.60.Cd 



I. INTRODUCTION 

Particle oscillations are fundamental for a wide range of in- 
teresting physics: quark mixing by the Cabbibo-Kobayashi- 
Maskawa matrix, its leptonic analog for massive neutrinos 
lH hypothetical photon-axion and photon-graviton oscil- 
lations in the presence of external magnetic fields ||3l|4|], and 
K° - K° oscillations |5]. Physically, these quantum sys- 
tems are coupled to macroscopic systems and through ex- 
ternal interaction their quantum evolution is altered. One 
prominent astrophysical context in which such oscillation for 
macroscopic systems is important involves the neutrinos in 
supernova cores that may execute flavor oscillations while 
simultaneously interacting with ambient supernova matter 
llal7ll8ll9l lirllllll . The primary motivation of the present pa- 
per is to provide a straightforward generalization of the Boltz- 
mann formalism with which to analyze the kinetics of oscil- 
lating neutrinos with collisions. By taking ensemble-averaged 
matrix elements of quantum field operators for mixed parti- 
cles, following the pioneering work of | I^ fl3lll^fl5lflql . we 
obtain quasi-classical phase-space densities that satisfy real- 
valued Boltzmann equations with coupling terms that account 
for the neutrino oscillations. The formalism is clear, numeri- 
cally tractable, and does not contain operators, complex quan- 
tities, or wavefunctions. 

In |nj we introduce the Wigner phase-space density oper- 
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ator approach from which we derive our formalism involv- 
ing classical phase-space neutrino flavor densities and their 
off-diagonal, overlap correlates. The latter couple the dif- 
ferent flavor states to account for neutrino oscillations. In 
SIIIII we demonstrate with several simple examples that the 
set of equations reproduces (i) standard flavor oscillations in 
a vacuum, (ii) flavor oscillations with absorptive matter cou- 
pling ("quantum decoherence" fl7io . and (iii) resonant matter- 
induced flavor conversion (the MSW effect) for a neutrino 
beam fisl Il9l Eflll . In flIVI we summarize the salient fea- 
tures of our new approach. We show in AppendixlAlhow one 
might include neutrino and antineutrino self-interactions into 
our formalism. 



II. GENERALIZED BOLTZMANN EQUATIONS 

To analyze a multi-particle system and reflect its inherent 
statistics, one usually sets up the density matrix of the system 
For classical systems this is the phase-space density. 
For quantum fields, the conceptually most similar analog is 
the Wigner phase-space density operator 1I2H.I22I1 

p(r, p, t) = f -^- e -* r yt( r - I r ', r)<A(r + \r' , t) , 
J (27T«) 2 2 

(1) 

where ip* and ifr denote the creation and annihilation opera- 
tors, respectively, for the mixing particles of interest lll2lll6ll . 
p works for fermions or bosons (resulting in different commu- 
tation relations of the creation and annihilation operators) and 
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massive or massless particles. Here, we average over an en- 
semble of two neutrino species, for example v e 's and v^'s and 
take the matrix elements in the number-density basis of Fock 
space: 



r = (nMnj) = (j: f r 



(2) 



where the indices i, j run over the neutrino flavor, and * 
means complex conjugation. A generalization to more par- 
ticle species is straightforward. The diagonal terms are real- 
valued and denote quasi-classical phase-space densities. The 
off-diagonal entries are complex-valued macroscopic over- 
lap functions. For completely decohered ensembles and for 
non-mixing ensembles the off-diagonal entries vanish. A 
Heisenberg-Boltzmann-type equation P fl^. Fnl El I2H1 
which can heuristically be derived by merging the Boltzmann 
equation with the Heisenberg equation, governs the time evo- 
lution of the matrix elements of the Wigner phase-space den- 
sity operator: 



dT 1 
~~dt + 2 



v, 



dT 
dr 



1 



dT 
dp 



= -i[Cl,T]+C , 



(3) 



where curly brackets denote matrix anti-commutators, square- 
brackets are matrix commutators, C is the collision matrix and 
Q is the mixing Hamiltonian. We componentize this equation 
for two oscillating neutrino species interacting with a back- 
ground medium. Then, 



C = 



C Ve 
C v , 



(4) 



and the following general mixing Hamiltonian reads: 
D(r, p, t) = Q(e, r) + Q vv (r, p, f) - Q V y(r, p, f) , 



(5) 



where Q(s, r) encompasses the vacuum mixing and the 
neutrino-matter interaction amplitude for the v e neutrino 
connected with the matter-induced mass. Q vy is the off- 
diagonal mixing contribution from neutrino-neutrino self- 
interactions and O vv is the analog for neutrino-antineutrino in- 
teractions. We show how to include neutrino self-interactions 
and neutrino-antineutrino interaction into our formalism in 
Appendix |X] and work here with the vacuum and ordinary 
matter contribution only. Therefore, we have for Q(e, r) the 
expression 



D(e, r) = 



nc - cos 26 + 2A sin 26 



L \ sin 26 cos 26 



(6) 



where 6 is the neutrino vacuum mixing angle between v e and 
and A is given by lEMEfill . 



. (L\2^G F 

A = — — - n e {r) 

\jtc I n 



(7) 



where L is the vacuum neutrino oscillation length: 

Anhce 



Am 2 c 



(8) 



Gf denotes Fermi's constant, n e (r) the electron number den- 
sity, s the neutrino energy, ni\ and mi are the masses of the 
neutrino mass eigenstates, and Am 2 = - m 2 . The other 
variables have their standard meanings. For antineutrinos the 
sign of A is reversed. Note that A, through n e , is in general 
a function of spatial position. In this paper, we do not take 
into account the effects of microscopic density fluctuations on 
the mixing in the ensemble l27ll . In other words, we assume 
that the scale of the spatial variation of the phase-space den- 
sity is larger than the neutrino de Broglie wavelength. Simi- 
larly, we assume small external forces. With these reasonable 
assumptions, we can ignore the off-diagonal terms on the left- 
hand-side of Eq. 0. Defining the real part of the off-diagonal 
macroscopic overlap in Eq. as 



fr — 2 (f e f + fefl) ' 

and the corresponding imaginary part as 



(9) 



(10) 



respectively, we find the generalized Boltzmann equations for 
real-valued quantities: 



df Ve df Ve . dfy^ 
dt +V ' dr +P ' dp 



2nc 

— —fiSm26 + C Ve 



dfr df r . df r 

dt dr F dp 



L 

2nc 



ft (cos 26 - A) 



dfi ... d fi _ Bfi _ 2nc^f Ve -U sin2g 



dt + V ' dr + P ' dp 



L v 2 
f r (cos 26 - A) ) . 



(11) 



The neutrino oscillations are incorporated with new sink 
and source terms indirectly coupling the standard Boltzmann 
equations for f, e and f Y through the off-diagonal macroscopic 
overlap functions /,- and f. The number of equations is in- 
creased from two to four. In principle, one can insert blocking 
factors for both the collision terms and the oscillation terms. 
Note that in the absence of collisions, Liouville's theorem for 
the total lepton specific intensity, dldt{j Vt +/ V J = 0, is re- 
covered. 

It is instructive to rewrite Eq. dl 11 1 in terms of the specific 
intensities of the neutrino radiation field. We use the one-to- 
one relation [28, 29] between the invariant phase-space den- 
sities, f Ve , f Vfl and the specific intensities, I Ve , I v , and define 
accordingly: 



L = 



(2nn) 3 c 2 



(2nK) 3 c 2 



(12) 
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For the off-diagonal macroscopic overlap analogs one defines 



(2nh) 3 c 2 



J, 



en 



(2ttH) 3 c 2 



(13) 



J 



The generalized Boltzmann equations in the laboratory (Eule- 
rian frame) for the radiation field of two oscillating neutrino 
species with collisions are then: 



\ dI Vt v dI Ve e 3 p <5(/ Ve e 3 ) 

c dt c dr c dp 

c dt c dr c dp 

\&R eil | y &R eil | g 3 p d(ft eti s-i) 

c dt c dr c dp 



2n 

-—I efl sm26 + C Ve 



2n . 



-I eu sin20 + C v 



~(cos26-A)I eft 



\dl eil v dl elx e 3 p d{l ell e 3 ) 2n(I Ve -L 



c dt c dr 



dp 



sin26» + (cos26»-A)K 



(14) 



r 



These equations and Eq. il Q are our major results. The colli- 
sion terms can be conveniently written as [28] 



C' 



-K S I + K U 



j <& Ve (n,n')iv e (&')d& , (15) 

where cD Vf is a phase function for scattering into the beam in- 
tegrated over the solid angle d£l', K% is the sum of all absorp- 
tion processes 2,- «;0"", where n,- is the number density of mat- 
ter species i and <r° denotes the absorption cross sections (for 
scattering processes the superscript a is replaced with s), TyJ 
is the equilibrium Fermi-Dirac occupation probability, and B Ve 
is the corresponding blackbody specific intensity. Changing 
the subscript from v e to yields the corresponding parame- 
ters for the v^'s. For sterile neutrinos, one substitutes v, for v M 
and sets the scattering and absorption terms to zero. Similarly, 
one can write a set of equations for antineutrinos with differ- 
ent collision terms and with the sign of A reversed. Neutrino 
and antineutrino evolution are implicitly coupled through pair 
processes. If self-interactions are included as we show in Ap- 
pendix|X] neutrino and antineutrino evolution are nonlinearly 
coupled. 



III. SIMPLE TESTS OF THE NEW FORMALISM 



oscillations are suppressed UHl . To demonstrate the expected 
limiting behavior of our formalism we artificially "turn off" 
matter suppression by setting the matter term A to zero. This 
is not the situation found in nature. We define the approximate 
oscillation time 



tnsc — 



L _ The 
2nc Am 2 c 4 



(16) 



and set the interaction rate of the v e 's equal to the interaction 
rate of the v^'s to define the characteristic absorption time 



t , = 

c " 1 Cc 



cN A pY n crf in 



(17) 



Na denotes Avogadro's number, Y„ is the neutron fraction per 
nucleon, and cr" v n is the cross section for absorption on neu- 
trons (see AppendixlEV We define the ratio of the oscillation 
to the absorption timescale: 



a = 



(18) 



col 



A. Oscillations with absorptive matter coupling 



and the dimensionless time coordinate t: 



It is straightforward to show that our set of generalized cou- 
pled Boltzmann equations behaves as we would expect from 
our experience with the standard wave function approach. As 
our first example, we solve Eq. (I14> for v e - v M oscillations in 
box of isotropic neutrinos that can also experience decohering 
absorption on matter. Note that in reality neutrino interac- 
tions only play a role at densities where observable neutrino 



(19) 



Furthermore, we set B v = B v and denote the dimensionless 
specific intensities and the off-diagonal macroscopic overlap 
functions that are normalized to the blackbody intensity with 
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—ley. Sm 2® ' 



a hat. The resulting dimensionless version of Eq. dl4> reads: 

d% 
8t 



sin20 + ar(l - 7 V(1 ) 



<3t 

~~dT 



—leu COS ^ 



sin 20 



■ Heii COS 20 



(20) 



where the scattering in and out of the beam has been canceled 
due to the assumption of isotropy. In a vacuum (a - 0) and 
with the initial conditions 7 Ve = 1, I v = 0, and consequently 
no off-diagonal overlap terms at r — 0, we obtain 



ZS.T) = 


1 - sin 2 20 sin 2 




%,{T) = 


sin 2 26 sin 2 f 




%(T) = 


- sin 20 cos 20 (cos t - 1 ) 
2 




I eli (T) = 


- sin 20 sin r . 

2 


(21) 



This behavior of the radiation field is unambiguously iden- 
tical to the probability density obtained by squaring the am- 
plitude of a single-neutrino wave function in a beam. The 
off-diagonal terms representing the macroscopic overlap peak 
when mixing of v e and v M neutrinos is maximal and vanish 
when the ensemble is single-flavored. In matter and for the 
initial conditions I v = 9i eM = I efl - and I Ve + 0, one can 
derive an harmonic oscillator equation for the early rate of 
evolution of 7 V : 



d% 

8t 2 



+ I — - a L I I Ve — const . 



(22) 



As expected, for a <K 1, the early time dependence of the so- 
lution is predominantly sinusoidal: neutrino oscillations dom- 
inate. For a » 1, an exponential decay/increase dominates. 
The timescale then is I /a. 

In Fig. ^ we depict the solutions to Eq. d20l for os- 
cillating v e 's in a box with nucleons. Initial conditions are 
I y = 'Reji = I en = and I Ye = 0.8 . For this example, flavor 
oscillations and collisions happen on the same timescale. The 
ensemble is guided to flavor and radiative equilibrium. Coher- 
ent flavor oscillations are disrupted by absorptive collisions. 
Asymptotically, the diagonal specific intensities for the v e 's 
and v^'s equilibrate. Absorption on neutrons, and by detailed 
balance, the resulting emissivity, drive the v e 's to the black- 
body intensity. The oscillation amplitude decreases with time; 
the quantum evolution of the system is decohered through ab- 
sorptive coupling with matter ("quantum decoherence" I17I0 . 
The real part of the off-diagonal overlap, ( R e ^ takes predomi- 
nantly negative values whereas the imaginary part, I~ CjU , oscil- 
lates symmetrically around zero. Both vanish asymptotically 
and no oscillations persist. 
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FIG. 1 : v e - oscillations of isotropic neutrinos in a box with ab- 
sorptive matter coupling, (a) Specific intensities, (b) Off-diagonal 
macroscopic overlap functions. Parameters: s Ve = s v = 10 MeV, 
p = 8 x 10 12 g citT 3 , T = 5 MeV, and from the large-mixing-angle 
solution (LMA)|3C|l: sin 26 = 0.9, and Am 2 = 6.9 x 10~ 5 eV 2 . A is 
artificially set to zero. 



B. Matter-enhanced resonant flavor conversion 



To demonstrate that our formalism contains the MSW ef- 
fect fljl fl9l Eoll . we solve Eq. H41 for a monoenergetic 
one-dimensional neutrino beam propagating down a density 
profile for which resonant matter-enhanced flavor conversion 
takes place. We define the dimensionless distance coordinate 
in terms of the oscillation length: 



2n 



(23) 



and the dimensionless matter-induced mass term in terms of 
its resonance value: 



A = 



A res cos 20 



(24) 
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FIG. 2: MSW effect, (a) Specific intensities, 
macroscopic overlap functions. An initially v e 
down the density profile &, where /? = 900. 



(b) Off-diagonal 
beam propagates 



The beam passes the resonance density for A = 1. We analyze 
the following dimensionless version of Eq. (I14> : 



dI Ve 
dx 

dx 

dx 

dx 



-I efl sin 26 



I eil sin 26 



= -cos26»(l -A) I m 



— sin 26 + cos 29 (l -A)<R ei i 



(25) 



where C' and C' have been set to zero. 

In Fig. |2J we depict the solutions to Eq. J25l > for a density 
profile of A = A, We set ft = 900 and thus ensure that the 
scale of spatial inhomogeneities is large compared to the mi- 
croscopic length scales such as the neutrino de Broglie wave- 
length and the oscillation length. Initially, the beam contains 
only v e neutrinos. The mixing angle is arbitrarily taken to be 
sin 2 26 = 0.18 (at present the LMA is favored 1 30]). From Fig. 
12 it is clear that at the resonance density, jc = -yjS = 30, the 
flavor composition of the beam is radically altered. For higher 
values of % the beam executes vacuum oscillations. In this 



illustrative problem the density at production is much greater 
than the resonance density. Then, the spatially averaged sur- 
vival probability of a v e neutrino going from matter to free 
space should be jH 



(P(v e _» v ,)> « (1 - P x ) sin 2 6 + P x cos 2 



(26) 



where P x is the Landau-Zener probability for non-adiabatic 
transitions. For the chosen density profile, propagation is adi- 
abatic and P x = 0. The averaged v e survival probability in 
Fig. |3]converges toward ~ 0.05, which is congruent with the 
value predicted in Eq. \26V . sin 2 6 ~ 0.05. For high densities 
(x < 20), matter suppression is severe. No flavor oscillations 
happen and the off-diagonal overlap functions 7?^ and I ' ejl are 
close to zero. In free space for ~x > 60, the imaginary part I ei , 
oscillates symmetrically around zero whereas the real part 'R^ 
is positive. 

The solution given in Fig. (0 is numerically equivalent to 
the solution obtained using the standard wavefunction formal- 
ism ElfilEolEi: 



(27) 



l dx\i/fyJ~2\ sin26» cos26»j\^vJ 



when we identify the specific intensities with the probability 
densities 



K <-» l«AvJ 2 
h ~ l<Av„l 2 > 

and the macroscopic overlap functions with 

Ren « +<P* v ^v) 



(28) 



(29) 



Thus, our Boltzmann formalism is completely consistent with 
the existing description. 



IV. CONCLUSIONS 

In this paper, we have derived a generalized set of Boltz- 
mann equations for real-valued phase-space densities of os- 
cillating neutrinos interacting with a background medium. 
The off-diagonal functions of the Wigner phase-space den- 
sity matrix representing macroscopic overlap are explicitly 
included and serve to couple the flavor states to reflect neu- 
trino oscillation physics. Conceptually, we have reduced the 
time evolution of creation and annihilation operators to that 
of real valued phase-space densities without losing quantum- 
physical accuracy. Important quantum effects such as matter- 
enhanced resonant flavor conversion and "decoherence" fl7ll 
through matter coupling are correctly incorporated. The gen- 
eralized Boltzmann equations are simple and very similar to 
the equations of classical transport theory. Neutrino oscilla- 
tions are incorporated by new sink and source terms that indi- 
rectly couple the expanded set of equations. We have shown 



6 



how to include neutrino self-interactions in our formalism. 
The self-interactions nonlinearly couple neutrino and antineu- 
trino evolution and therefore for the most generic case one 
has to deal with eight non-trivially coupled equations for a 
two-flavor ensemble and their antiparticles. 

Using this formalism, codes that now solve the standard 
Boltzmann equations for the classical neutrino phase-space 
density (/ Vf ), or which address its angular and/or energy mo- 
ments, can straightforwardly be reconfigured by the simple 
addition of source terms and similar transport equations for 
overlap densities that have the same units as f V) , to incorpo- 
rate neutrino oscillations in a quantum-physic ally consistent 
fashion. 
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trino with momentum q and the coupling strength is given by 

B = ^ Gf . The neutrino fields of the background neutrino 
with momentum q are normalized such that 



I+ K) 



1 



(A3) 



We now match the above expressions and convert them into 
our quasi-classical phase-space densities and the macroscopic 
overlap functions. The v-v mixing Hamiltonian for test neutri- 
nos with momentum p in ensemble-averaged form is denoted 
by 



<-» i B Vr B r + iB. 

Q w (p,r,0 = U_ k . 



(A4) 



where, while throwing away the overall phase term in Eq. 
(IA2> proportional to the identity matrix, we have for the di- 
agonal elements: 

B Vt (p,r,0 = B j d 3 q(l - cos6» pq )/ Vf (q,r,r) 
Bv„(p,r,f) = pf ^q(l-cosfl")/ v ,,(q,i\f). 

(A5) 



APPENDIX A: SELF-INTERACTIONS 

Evolution in a neutrino background is nontrivial l32l l33l 
l34ll35ll . To address this, one can retain the one-particle char- 
acter of description and neglect any effects arising from co- 
herent many^ody state formation or flavor entanglement l36l 
y7|, |3]||39|, 40]. The local low-energy four Fermi-interaction 
Hamiltonian of neutrino-neutrino scattering is given by 



Q 1 



pq 



2l 

V2 



(Al) 



where the sum is over all neutrino flavors. This effective 
Hamiltonian must satisfy U(N) flavor symmetry for a system 
of N flavors |33f. For a two-flavor system consisting of v c 's 
and v^'s one can rewrite this Hamiltonian to accentuate its off- 



I , q I 2 , q , q*^ 

Wv,\ <PvAv u 



diagonal character 11341 13 

Q p v q = B{\ -cos<9 pq )x 



(A2) 



where the coupling coefficient includes the angle between the 
"test neutrino" with momentum p and the background neu- 



We used the "matching" from our formalism to the wavefunc- 
tion formalism as prescribed in Eqs. d28l > and ( 12 9> . The mo- 
mentum integration goes over all momenta q in the ensemble. 
For the off-diagonal elements we write in our notation, 

fi,-(p,r,f) = B j d\(l -cosflW)/ r (q,r,r) 
Bi(p,r,0 = B J d 3 q(l - cos fl") /-(q, r, f) , 

(A6) 

where we have used the variables of Eqs. (|9jl and dlOi . For the 
neutrino-antineutrino interaction there is 



Q vv (p, r, f) 



B r - iBi 



B r + iBi 

B v ., 



(Al) 



where the B's are now defined in terms of the antineutrino 
phase-space densities f Ve , f v , f r , and Note that in accor- 
dance with Eq. Q the coupling coefficient in front of the 
integrals has to be implemented in the equations for neutri- 
nos with reversed sign. For two mixing neutrino species and 
their antiparticles interacting with a background medium and 
with neutrino-neutrino interactions included, the generalized 
Boltzmann equations in their most generic form are: 
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dfv e df^ . dfy^ 
dt V ' dr P dp 



df v , df. 



+ V- — + p 

dr dr F dp 



-fi 

fi 



^ sin 26 + 2 (B r + B r ) + 2f r (fi, + B,) + C v 



dt dr F dp 



L 
2nc 
~L 

^ (A - cos 26») + B Ve - B Ve - B Vfi + B v 



sin 26 + 2 (B r + B r ) - 2f r (fi, + £,) + C V(J 



(/v e -/v,)[ysin20 + (B r - J B r ) 



+ (/v.-/vJ(«i-B,) 



(A - cos 26) + B Ye - B Ye - By + By 



(A8) 



In a very similar fashion one can write the corresponding anti-particle analog. We need to interchange tildes and reverse the sign 
for A to complete our set of equations: 



dfv, df Ve . df Ve ~ 
— + v — + p- — = -fi 
dt dr dp 

^r +v -^r +p -^p- = fi 



2nc 



sin 20 + 2 (fi r + B r ) + 2f r (fi; + B,-) + C Ve 



2nc 



— sin 26 + 2 (b r + B r ) - 2f r (5,- + B,) + C v „ 



•9/r 5/ r . 3/ r - 

M +V .M + P .M 

ft dr F 5p 



27TC 

— ("A - COS 20) + By e - By e - By^ + By^ 



= (/v e -/v„)[ysin20 + (B r -B r ) 



— (-A - cos 26) + B Ve - B Ve - B Vfi + B v 



(A9) 



Neutrino and antineutrino evolution are nonlinearly coupled. The collision terms differ for neutrinos and antineutrinos. The sign 
reversal of A means that, dependent on the mass hierarchy, only neutrinos or antineutrinos execute the MSW resonance. The 
above nontrivially coupled set of eight equations is entirely real-valued; yet they contain all the quantum-mechanical oscillation 
phenomenology. This is the complete set of kinetic equations that include neutrino self-interactions for the real- valued neutrino 
phase-space densities f Ye and f v and the corresponding antineutrino phase-space densities f Vt and f v . 



by 



APPENDIX B: CROSS SECTION: v e + n — > e + p 

A convenient reference neutrino cross section is cr , given 



4G 2 F (m e c 2 ) 2 

7l(fic) 4 



1.705 x 10 _44 cm 2 , 



(Bl) 



the reaction v e + n — > e + p is then given by 



(l+3g 2 A(s Ve +A 



l - 



1/2 



By, 



(B2) 



where Gf is the Fermi weak coupling constant 1.436 x where gA is the axial-vector coupling constant ( 1.23), and 

10 -49 ergs cm 3 ). The total v e - n absorption cross section for A np = m n c 2 - m p c 2 = 1.29332 MeV . 
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